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摘要 边缘 检测 是 图 像 处 理 和 计算 机 视觉 中 的 基本 问题 ， 其 目的 是 标识 数字 图 像 中 亮度 变化 明显 的 点 。 激 
光 测 距 是 测 距 领域 精度 最 高 的 的 重要 技术 手段 ， 激 光 测 距 数据 能 够 被 应 用 于 科学 研究 之 前 ， 需 要 对 原始 
数据 进行 一 系列 的 预 处 理 。 根 据 目 标 信 号 和 噪声 信号 在 图 像 上 的 分 布 特点 ， 本 文 提出 了 一 种 基于 边缘 检 
测 来 提取 目标 信号 的 方法 ， 将 原始 数据 转化 成 测 距 残 差 图 像 ， 对 其 进行 边缘 检测 后 识别 信号 位 置 。 用 过 
对 实测 数据 处 理 表 明 本 文 提出 的 算法 具有 一 定 的 可 行 性 、 通 用 性 。 
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Abstract Edge detection is a basic problem in image processing and computer vision. Its purpose is to identify 
points with obvious brightness changes in digital images. Laser ranging is an important technical method with the 
highest accuracy in the field of ranging. Before laser ranging data can be used in scientific research, a series of 
preprocessing of the original data is required. According to the distribution characteristics of the target signal and 
noise signal on the image, this paper proposes a method to extract the target signal based on edge detection, which 
converts the original data into the residual image of distance measurement, and then recognizes the signal position 
after edge detection. It has been used to process the measured data that the algorithm proposed in this paper has a 
certain feasibility and versatility. 
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激光 测 距 是 当前 获取 目标 高 精度 距离 的 的 重要 技术 手段 之 一 ， 在 监测 大 陆 板 块 运动 、 地 过 形变 、 地 
球 上 自转 、 改 进 地球 重 力 场 和 地 心 引 力 参数 ， 确 定 地 球 和 海洋 潮汐 变化 的 规律 ， 监 测 空间 碎片 等 方面 具有 
重要 作用 ”，。 随 着 商业 航天 近 些 年 的 大 力 发 展 ， 空 间 目标 的 数量 也 呈 和 急剧 增长 ， 我 们 更 加 需要 对 空间 目 
标 有 效 的 高 精度 测量 。 测 距 数 据 能 够 被 应 用 于 各 领域 之 前 ， 需 要 对 原始 数据 进行 预 处 理 。 

图 像 的 边缘 是 图 像 最 基本 也 是 最 重要 的 特征 之 一 ， 边 缘 检测 的 目的 是 去 发 现 图 像 中 关于 形状 和 反射 
或 透射 比 的 信息 “。 目 前 卫星 激光 测 距 信号 提取 的 方法 主要 有 人 工 识 别 ，Graz 回 波 自动 识别 ， 还 有 
poisson 滤波 算法 等 。 针 对 以 上 算法 需要 设置 过 多 参数 和 人 为 不 确定 性 因素 的 不 足 ， 本 文 提 出 了 一 种 基于 
图 像 边缘 检测 的 激光 测 距 信号 提取 方法 ， 并 结合 了 和 一 些 常规 的 数据 处 理 方法 ， 提 取出 卫星 有 效 回 波 信 


叮 


To 


2 基本 原理 


卫星 激光 测 距 测量 的 是 激光 脉冲 往返 的 飞行 时 间 趾 ， 在 测 距 过 程 中 ， 由 于 系统 噪声 、 环 境 噪声 等 的 影 
响 ， 残 差 数 据 中 有 大 量 的 噪声 ， 我 们 需要 从 噪声 中 识别 信号 数据 。 本 文 处 理 的 数据 对 象 是 测 距 残 差 图 
像 ， 我 们 对 残 差 图 像 进 行 卷 积 操作 ， 得 到 图 像 的 边缘 信息 ， 即 得 到 信号 位 置 ， 然 后 再 通过 图 像 的 对 应 关 
系 还 原 到 原始 信号 ， 这 一 阶段 的 数据 含有 比较 多 的 噪声 数据 ， 我 们 仍 需要 经 过 对 数据 进行 拟 合 去 噪 得 到 
需要 的 信号 数据 ， 根 据 其 体 精度 的 要 求 可 以 继续 对 数据 进行 处 理 。 


3 方法 步骤 


3.1 获 取 测 距 残 差 图 像 

本 次 原始 测 距 数据 全 部 源 自 于 云南 天 文 台 应 用 天 文 研究 团 组 1.2 米 望远镜 ， 观 测 目标 为 卫星 和 部 分 空 
间 碎 片 等 。 激 光 器 百 赫兹 或 其 他 不 同 频率 都 需要 插值 计算 激光 器 每 个 主 波 时 刻 对 应 的 预报 距离 趾 。 本 文采 
用 的 是 拉 格 朗 日 三 阶 插值 。 目 前 ， 卫 星 测 距 系统 一 般 都 采用 事件 计时 器 作为 时 间 测量 设备 。 该 器 件 将 每 
个 激光 脉冲 的 发 射 时 刻 《〈 主 波 ) 与 接收 时 刻 〈 回 波 ) 分 别 作为 A、B 事件 ， 需 要 将 回 波 与 真实 的 主 波 进行 
匹配 ， 匹 配 过 程 的 流程 是 ， 对 于 每 个 B 事件 (bj)， 与 每 一 个 A 事件 〈ai) 的 时 刻 做 减法 ， 其 差 值 再 减 去 
此 A 事件 所 对 应 的 预报 距离 时 间 ， 当 其 小 于 给 定 的 阔 值 T 时 则 认为 匹配 成 功 器。 


|(b;—ai)— pred(ai) k=7T (1) 


202101.00040v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


遍历 公式 (1) 中 的 ij， 满 足 条 件 的 保存 主 波 时 刻 和 所 匹配 成 功 的 回 波 时 刻 ， 然 后 生成 测 距 残 差 
像 ， 横 坐标 是 主 波 时 刻 Ta， 纵 坐标 是 残 差 值 。 
3.2 边缘 检测 

图 像 边 缘 具 有 方向 和 幅度 两 个 特征 ， 沿 着 图 像 边缘 方向 ， 像 素 的 亮度 变化 较 小 ， 但 在 垂直 于 边缘 方 
向 ， 像 素 的 亮度 值 变化 比较 剧烈 @@。 图 像 边缘 检测 的 基本 思路 就 是 通过 对 图 像 求解 一 阶 导 数 的 局 部 极 大 
值 来 表示 图 片 的 边缘 信息 。 经 典 的 边缘 检测 方法 是 对 原始 图 像 中 所 有 像素 的 周围 邻 域 进行 卷 积 运算 ， 这 
种 运算 的 模块 贝 称 为 检测 算 子 ， 和 常用 的 边缘 检测 算 子 有 : Roberts 算 子 、Sobel 算 子 、Prewitt 算 子 、Kirsch 
算 子 、LOG 〔〈 高 斯 一 拉 普 拉 斯 ) 算 子 和 Canny 算 子 等 。 这 些 方法 多 是 以 要 人 处理 像素 为 中 心 的 邻 域 进行 灰 
度 分 析 ， 实 现 对 图 像 边缘 的 提取 品 。 


碳 


下 面 介绍 常用 的 一 种 边缘 检测 模板 Sobel 算 子 。Sobel 算 子 主要 是 利用 像素 点 上 下 、 左 右 邻 点 的 灰 度 
加 权 算 法 ， 县 有 很 好 的 边缘 检测 性 ， 并 对 噪声 具有 平滑 作用 ， 缺 点 是 它 同时 也 会 检测 出 许多 的 伪 边 缘 。 
一 般 来 说 对 精度 要 求 不 是 很 高 时 ， 经 常 采用 这 种 方法 。Sobel 算 子 有 两 个 : 一 个 是 检测 水 平 边 缘 的 ， 另 一 
个 是 检测 垂直 边缘 的 。Sobel 算 子 是 滤波 算 子 的 形式 ， 用 于 提取 边缘 ， 可 以 利用 快速 卷 积 函 数 ， 简单 有 
效 ， 因 此 应 用 广泛 。Sobel 算 子 边缘 算 子 所 采用 的 算法 是 先进 行 加 权 平 均 ， 然 后 进行 微分 运算 ， 算 子 的 计 
算 方法 如 下 : 


Gx=[fCx-1, y-D+2fCx-l, y+f x-l 7+D]-[Fx+D y-D+2fGx+tl, +flxtl, y+D| 2) 


Gy=[f Cx-l, y+D+2fCx, y+D+f xtl, y+D]-[fCx-l y-D+2fCx y-D+f extl, >-D] 3) 


公式 (2) 是 检测 水 平 边 缘 ， 公 式 〈3) 是 检测 垂直 边缘 。 


算 子 包含 两 组 3x3 的 和 矩阵， 分 别 为 横向 及 纵向 ， 将 之 与 图 像 作 平 面 卷 积 ， 即 可 分 别 得 出 横向 及 纵向 
的 亮度 差分 近似 值 。 如 果 以 A 代表 原始 图 像 ，Gx 及 Gy 分 别 代表 经 横向 及 纵向 边缘 检测 的 图 像 ， 其 方法 
具体 如 下 ， 其 中 I 是 图 像 G 的 一 部 分 ， 循 环 遍历 图 像 G 得 到 经 过 边缘 检测 之 后 的 图 像 。 


水 平 变化 : 将 I 与 一 个 奇数 大 小 的 内 核 G, 进行 卷 积 。 当 内 核 大 小 为 3 时 , G; 的 计算 结果 如 公式 
(4) 。 


1 0 -| 
Or-|2 0 -21*7 (4) 
1 0 -1l 


垂直 变化 : 将 I 与 一 个 奇数 大 小 的 内 核 Gy 进行 卷 积 。 当 内 核 大 小 为 3 时 , G， 的 计算 结果 公式 (5) 。 


1 2 1 
G=I0 0 O04*7 ($5) 
-1 -2 -| 


我 们 对 已 有 残 差 图 像 进行 基于 Sobel 算 子 的 边缘 检测 ， 将 以 2019 年 11 月 15 日 glonass-129 卫星 为 


例 ， 如 图 1 和 图 2。 
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图 1 2019.11.15-glonass129 测 距 残 差 原 图 


Fig. 1 2019.11.15-glonass129 original image of ranging residual 


图 2 2019.11.15-glonass129 残 差 图 经 过 sobel 方法 边缘 检测 


Fig. 2 2019.11.15-glonass129 residual image after Sobel method edge detection 


为 了 方便 处 理 ， 我 们 将 原 图 转 成 灰 度 图 进行 sobel 边缘 检测 处 理 ， 结 果 如 图 2， 这 种 效果 不 是 我 们 想 


要 的 ， 从 图 1 中 可 以 看 到 信号 集中 在 中 间 的 黑色 宽 线 段 ， 然 而 经 过 sobel 边缘 检测 后 ， 图 1 中 所 有 的 位 置 


边缘 信息 全 部 检测 并 保留 下 来 ， 但 是 我 们 并 不 需要 那些 孤立 的 噪声 点 ， 因 此 ， 根 据 sobel 背后 的 数学 原 


AH 


， 我 们 设计 出 了 一 种 专门 针对 测 距 信号 的 边缘 检测 ， 也 就 是 我 们 需要 设计 一 个 大 尺度 的 卷 积 核 〈 算 
子 ) 并 行 近似 横向 的 识别 ， 因 为 在 大 斥 度 的 残 差 图 中 ， 信 和 号 数据 显示 在 残 差 图 上 一 般 是 近似 的 直 


线 。 这 些 参 数 的 选择 需要 尽 可 能 的 包含 信号 ， 当 然 会 把 噪声 也 涵 善 进来 ， 由 于 在 残 差 图 像 信号 位 置 附近 
信和 号 密度 一 般 比 噪声 密度 大 ， 后 面具 需 简 单 的 经 过 统计 学 的 原则 和 粗 拟 合 来 去 除 过 于 明显 的 噪声 。 


sobel 算 子 的 尺寸 是 3*3， 九 个 元 素 对 图 像 进行 遍历 卷 积 ， 经 过 图 2 的 观察 分 析 ， 对 于 激光 测 距 残 关 
图 来 说 ， 小 尺寸 的 卷 积 操作 会 保留 更 多 的 噪声 。 因 此 ， 本 文 设计 了 尺寸 为 5*65 的 算 子 ON 〈 具 体形 式 见 
公式 6) ,来 对 测 距 残 差 图 像 进行 卷 积 操作 ， 图 3 是 经 过 新 算 子 的 处 理 结果 。 


1 1 11 
0 0.…00 

Ov-|0 0 .… 0 0|/65 (6) 
0 0.… 00 
EE i | 


公式 (6) 中 的 矩阵 大 小 是 5*65， 需 要 补充 的 是 所 有 的 测 距 残 差 图像 原 始 尺寸 均 为 1610*2425*3( 宽 
* 长 * 通 道 ) ， 最 后 边缘 检测 结果 的 图 片 尺 寸 均 为 1068*2361 〈 宽 * 长 ) ， 由 于 对 原始 图 片 进 行 了 灰 度 处 理 
所 以 结果 变 为 单 通道 。 下 面 说 明 参 数 大 小 设置 原因 ， 根 据 信号 长 度 的 特点 ， 经 过 对 一 些 图 片 处 理 之 后 ， 
设置 了 算 子 经 验 长 度 值 65; 结构 类 似 与 sobel 算 子 ， 第 一 行 与 第 三 行 是 相反 数 ， 第 二 行 设 置 为 零 ， 与 
sobel 不 同 的 是 新 的 算 子 第 一 行 全 部 设置 为 1， 因为 在 图 片 转化 成 灰 度 图 后 ， 为 加 快速 度 方 便 计 算 ， 本 文 


又 将 其 二 值 化 ， 所 以 横向 的 像素 点 强度 服从 均匀 分 布 而 不 是 高 斯 分 布 ， 所 以 全 部 设置 为 一 样 ， 同 样 也 对 
和 矩阵 进行 了 零 均 值 和 归 一 化 处 理 。 
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图 3 2019.11.15-glonass129 残 差 图 经 过 超大 尺度 算 子 边缘 检测 结果 
Fig.3 2019.11.15-glonass129 residual image after super-large-scale operator edge detection results 


3.3 信和 号 恢复 


经 过 对 图 像 进行 边缘 检测 ， 
出 的 是 主流 和 回 波 已 经 在 第 一 
刻 ， 本 文采 用 的 做 法 是 按照 图 像 大 小 比例 和 信 


最 终 得 到 了 信号 在 图 像 
步 匹配 完成 ， 所 以 此 时 得 到 主 波 信号 


号 位 置 比 例 来 还 原 时 间 信号 


Ph 的 位 置 ， 
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但 是 我 们 需要 的 是 主 波 时 刻 ， 需 要 指 
即 可 。 所 以 对 图 像 处 理 还 原 到 主 波 时 


。 为 便于 描述 ， 将 原始 残 差 图 
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像 记 为 Go， 得 到 的 图 像 记 为 G1， 确定 Go 和 Gi 大 小 比例 系数 k，G 图 起 始 主 波 时 刻 Ao ， 有 具体 做 法 是 ， 
对 于 Gi 中 的 每 个 信号 点 的 横 化 标 Ci 所 对 应 的 原始 主 波 时 刻 计算 公式 为 : 


Ci= Amin|(A0+ Bi*hk)| (7) 


这 时 得 到 信号 为 粗略 信号 ， 这 意味 着 还 原 主 波 数据 中 还 会 携带 一 些 噪声 数据 如 图 4， 我 们 仍 需 要 过 滤 掉 明 


显 的 噪声 点 。 
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图 4 2019.11.15-glonass129 数据 还 原 后 含有 噪声 结果 


Fig. 4 2019.11.15-glonass129 data contains noisy results after restoration 
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图 5 glonass129 数据 还 原 后 粗略 去 除 噪声 结果 
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Fig. 5 glonass129 Data reduction results after rough reduction 


4.1 处 理 结果 


我 们 对 不 同 的 卫星 和 一 些 空间 碎片 都 进行 了 边缘 检测 ， 但 是 由 
标 来 说 明 ， 这 些 原始 信号 包括 弱 噪 声 强 信号 、 强 噪声 弱 信 号 、 强 噪声 强 信 号 。 如 
原始 残 差 图 ， 处 理 结果 (b) 是 边缘 检测 后 


Original Image 


Processed Image 


的 信号 残 差 图 ，(c) 放 大 显示 的 残 差 图 


于 篇 幅 有 限 ， 本 文 只 
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只 列举 一 些 典 型 
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型 目 
图 6-- 图 9， 其 中 (a) 图 是 
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图 6 37731 空间 碎片 的 原始 数据 (a), 处 理 结果 (b), (c) 
Fig. 6 37731 space debris raw data and processing results 
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图 7 beaconc 卫星 的 原始 数据 (a), 处 理 结果 (b), (c) 
Fig.7 beaconc satellite raw data and processing results 
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图 8 021108 空间 碎片 的 原始 数据 (a), 处 到 
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结果 (b), (c) 


Fig. 8 021108 Space debris raw data and processing results 
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图 9 09904 空间 碎片 的 原始 数据 (a), 处 理 结果 (b), (c) 
Fig. 9 09904 Space debris raw data and processing results 
表 1 一些 目标 的 精度 分 析 
Table 1 Precision analysis of some targets 
NORAD ID RCS/m” Precision/cm Number 
Edge_M Graz Edge_M Graz 
37731 8.5270 213.1 256.3 164 131 
021108 5.6696 285.2 330.9 94 52 
09904 7.7161 174.5 175.6 79 46 
1328(beaconc) 2.2651 40.4 39.7 1454 1417 
17908(ajisai) 3.9811 80.4 82.0 34021 35364 


表 1 是 一 些 空间 目标 的 精度 分 析 ， 第 一 列 NORAD ID 是 北美 防空 司令 部 空间 物体 的 编号 ， 第 二 列 是 
雷达 散射 截面 积 (Radar Cross Section，RCS)， 第 三 列 是 双 问 测 距 精度 (激光 器 脉 宽 7ns)， 子 列 Edge_M 
是 边缘 检测 方法 的 精度 ， 子 列 Graz 是 目前 常用 测 距 数据 处 理 方法 (Graz 自动 识别 方法 ) 的 精度 ， 第 四 列 是 


识别 信号 点 数量 ， 子 列 分 别 是 两 种 方法 识别 的 点 数 。 
4 2 结果 分 析 


通过 4.1 结果 可 以 看 出 ， 通 过 边缘 检测 的 方法 识别 的 信号 长 度 较 原 图 短 一 些 即 初始 识别 存在 漏 点 ， 主 
要 原因 是 卷 积 核算 阵 的 长 度 过 长 ， 但 是 这 是 不 可 避免 的 ， 检 测 大 尺度 横向 的 边缘 ， 只 能 采取 大 尺度 卷 积 
核 。 为 了 使 数据 量 不 减少 ， 采 取 了 一 种 简单 的 扩充 方式 ， 其 体 是 把 边缘 检测 每 段 结果 的 长 度 左右 扩大 一 
些 ， 只 要 其 残 差 符合 给 定 的 阔 值 即 亦 认为 是 信号 ， 后 期 需要 再 进行 拟 合 分 析 。 通过 表 1 可 以 看 出 ， 碎 片 
的 精度 大 致 在 米 级 ，beaconc 和 ajisai 激光 测 距 卫 星 的 精度 为 分 米 级 。 


通过 以 上 分 析 可 以 知道 ， 本 文 在 信号 数量 较 少 时 ， 并 不 能 检测 到 边缘 ， 也 就 是 该 方法 有 识别 最 少 信 
号 量 的 限制 ， 需 要 手动 调整 卷 积 核 的 大 小 以 适应 ， 这 是 后 续 努 力 的 方向 ， 另 外 ， 新 的 卷 积 核 计算 比 较 耗 
缆 时 间 也 是 接 下 来 改进 的 方向 。 


本 


本 文 提出 了 一 种 基于 通过 图 像 边缘 检测 识别 激光 测 距 信号 的 新 方法 ， 该 方法 处 理 的 数据 对 象 是 测 距 
残 差 图 像 ， 对 其 进行 图 像 边缘 检测 ， 利 用 新 的 大 尺度 卷 积 核 来 完成 特定 的 检测 目标 ， 最 后 根据 图 像 对 应 
数学 关系 还 原 回 主 波 和 回 波 信号 。 为 了 特殊 目标 的 需求 ， 自 己 设计 专用 的 卷 积 核 〈 算 子 )， 男 由 于 图 像 还 
原 时 误差 较 大 来 了 额外 的 噪声 ， 通 过 对 所 列 目标 的 精度 分 析 表 明 该 方法 的 有 效 性 。 
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